BATTERY CHARGE INDICATOR SUCH AS FOR AN 
IMPLANTABLE MEDICAL DEVICE 

TECHNICAL FIELD 

5 This patent application pertains generally to batteries, and more particularly, 

but not by way of limitation, to a battery charge indicator such as for an implantable 
medical device. 

BACKGROUND 

10 Implantable medical devices include, among other things, cardiac rhythm 

management (CRM) devices such as pacers, cardioverters, defibrillators, cardiac 
^synchronization therapy (CRT) devices, as well as combination devices that 
provide more than one of these therapy modalities to a subject. Such implantable 
devices are typically powered by a battery. When the battery's useful life has 

15 expired, the implanted device is typically explanted and replaced. Therefore, it is 
often useful to know how much battery energy has been used and/or how much 
battery energy remains. 

BRIEF DESCRIPTION OF THE DRAWINGS 
20 In the drawings, which are not necessarily drawn to scale, like numerals 

describe substantially similar components throughout the several views. Like 
numerals having different letter suffixes represent different instances of 
substantially similar components. The drawings illustrate generally, by way of 
example, but not by way of limitation, various embodiments discussed in the present 
25 document. 

FIG. 1 is a schematic diagram illustrating generally an example of portions 
of a system that includes a battery and that provides a battery capacity indication. 

FIG. 2 is a flow chart illustrating generally a technique for obtaining a 
reading of how much battery capacity has been used (or remains). 
30 FIG. 3 is a flow chart illustrating generally an alternative technique for 

obtaining a reading of how much battery capacity has been used (or remains). 
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FIG. 4 is a graph of average battery terminal voltage (in millivolts) vs. 
battery capacity (in milliampere-hours) for each of nine different substantially 
constant current loads. 

FIG. 5 is a graph of battery terminal voltage (millivolts) vs. load resistance 
(kilo-ohms), as derived from the data of FIG. 4 at various levels of battery capacity. 

FIG. 6 is a graph of battery terminal voltage (millivolts) vs. load resistance 
(kilo-ohms), in which the load resistance is plotted on a logarithmic scale. 

FIG. 7 is a battery terminal voltage (millivolts) vs. load impedance (kilo 
ohms) graph of four representative fitted curves. 

FIG. 8 is a graph of battery terminal voltage (millivolts) vs. capacity 
(milliampere-hours) using both predicted and actual data for a 32.4 kQ load 
impedance. 

FIG. 9 is a graph of the residuals (millivolts) vs. capacity (milliampere- 
hours) for the 32.4 kQ load predicted and actual data of FIG. 8. 

FIG. 10 is a graph of battery terminal voltage (millivolts) vs. capacity 
(milliampere-hours) for predicted and actual data for several higher resistance loads. 

FIG. 11 is a graph of the residuals (millivolts) vs. capacity (milliampere- 
hours) for the 100 kQ load predicted data (which was predicted using both 32.4 kQ 
load actual data and 64.9kQ load actual data) and actual data. 

FIG. 12 is a graph of capacity (milliampere-hours) vs. battery terminal 
voltage (volts), illustrating both the resulting cubic regression line and the lines 
indicating the 95% prediction interval. 

DETAILED DESCRIPTION 
In the following detailed description, reference is made to the accompanying 
drawings which form a part hereof, and in which is shown by way of illustration 
specific embodiments in which the invention may be practiced. These 
embodiments, which are also referred to herein as "examples," are described in 
sufficient detail to enable those skilled in the art to practice the invention, and it is to 
be understood that the embodiments may be combined, or that other embodiments 
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may be utilized and that structural, logical and electrical changes may be made 
without departing from the scope of the present invention. The following detailed 
description is, therefore, not to be taken in a limiting sense, and the scope of the 
present invention is defined by the appended claims and their equivalents. 
5 In this document, the terms "a" or "an" are used, as is common in patent 

documents, to include one or more than one. In this document, the term "or" is used 
to refer to a nonexclusive or, unless otherwise indicated. Furthermore, all 
publications, patents, and patent documents referred to in this document are 
incorporated by reference herein in their entirety, as though individually 
1 0 incorporated by reference. In the event of inconsistent usages between this 
documents and those documents so incorporated by reference, the usage in the 
incorporated reference(s) should be considered supplementary to that of this 
document; for irreconcilable inconsistencies, the usage in this document controls. 
FIG. 1 is a schematic diagram illustrating generally, by way of example, but 
1 5 not by way of limitation, one example of portions of a system 1 00. In this example, 
the system 100 includes a battery 102 for storing energy. The system 100 also 
provides a battery capacity indication. The battery capacity indication relates to 
how much of the stored energy has been used (such as by delivering an electrical 
current to an actual load impedance 104) or to how much of the stored energy 
20 remains within the battery 102. In one example, the actual load impedance 104 
represents the effective impedance presented by various circuits powered by the 
battery 102, such as a microprocessor and/or other various circuits in an implantable 
cardiac rhythm management (CRM) device. 

In FIG. 1, the system 100 includes a load impedance measuring circuit 106 
25 for measuring the impedance value of the actual load impedance 104. In this 

example, the load impedance measuring circuit 106 includes a voltage measuring 
circuit 108, to measure the terminal voltage across the battery 102, and a current 
measuring circuit 110, to measure the load current delivered to the actual load 
impedance 104 by the battery 102. 
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A microprocessor or similar processor circuit 112 is coupled to the load 
impedance measuring circuit 106, such as to receive the measured battery terminal 
voltage value and the measured load current value. The processor 112 calculates the 
measured load impedance value of the actual load impedance 104 by dividing the 
5 measured battery terminal voltage value by the measured load current value. Using 
the measured load impedance value as an input parameter to battery capacity 
characterizing data that was previously stored in a memory circuit 114, the 
processor circuit 112 provides a resulting battery capacity reading. In one example, 
the processor circuit 112 calculating the battery capacity reading is located in an 
10 implantable CRM device. In another example, the processor circuit 112 calculating 
the battery capacity reading is located in an external device that is communicatively 
couplable to the implantable CRM device, such as to receive the battery terminal 
voltage and/or load impedance readings and/or any other needed data upon which 
the battery capacity reading is based. 
15 In one example, the characterizing battery data stored in the memory circuit 

114 is obtained from a separate characterizing battery 116 (or several such batteries 
116) that was previously discharged by a characterizing load impedance 118 to 
obtain the characterizing battery capacity data. In one example, this includes using 
a battery manufacturer's published capacity vs. voltage data, resulting from the 
20 manufacturer's discharge of at least one characterizing battery 116. In other 
instances, such characterizing data is independently obtained, such as by 
discharging the characterizing battery 116 through the characterizing test load 
impedance 118 during a time period, measuring the terminal voltage of the 
characterizing battery 116 and the characterizing test load current at determinable 
25 times during the time period, and calculating the charge removed from the 
characterizing battery 116 during the time period, such as by exactly or 
approximately integrating the characterizing test load current. 

The separate at least one characterizing battery 116 is manufactured in the 
same manner as the battery 102, e.g., using substantially the same battery chemistry 
30 and substantially the same manufacturing process. However, the characterizing load 
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impedance 118 may differ in value from the actual load impedance 104. For 
example, to more quickly obtain the characterizing battery capacity data, the 
characterizing load impedance value may be substantially less than the actual load 
impedance value, which may take several years to discharge using the actual load 
impedance 104. Among other things, this document discusses how to translate 
between the value of the actual load impedance 104 and the characterizing battery 
capacity data (taken at a possibly different value of a characterizing load impedance 
118), such that a meaningful battery capacity reading can be made using the 
characterizing battery data. 

FIG. 2 is a flow chart illustrating generally, by way of example, but not by 
way of limitation, one technique for obtaining a reading of how much battery 
capacity has been used (or remains). At 200, the terminal voltage of the battery 102 
is measured, such as by using the voltage measurement circuit 108. At 202, a 
corresponding measurement of load current is measured, such as by using the 
current measurement circuit 110. At 204, an actual load impedance value is 
calculated, such as by dividing the measured battery terminal voltage by the 
measured load current using the processor 112. At 206, the measured battery 
terminal voltage is translated into an adjusted voltage, Vadj\ for use as an input 
parameter into previously stored battery characterizing data in the memory circuit 
114, such as capacity vs. voltage data obtained from a different characterizing 
battery 116 that is substantially identical in chemistry and manufacture to the battery 
102. In one example, the adjusted voltage, Vadj, is calculated according to Equation 
1. 



In Equation 1, R2 is the value of the characterizing test load impedance 118, Rl is 
the value of the actual load impedance 104, as calculated at 204, b is a coefficient 
obtained from the characterizing of the battery 116, and Vobs is the battery terminal 
voltage value of the battery 102, as measured at 200. 



i = Fofo + &*log( — 
Rl 




(1) 
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At 208, the adjusted voltage, Vadj, is used as an input parameter into 
previously stored battery capacity vs. voltage characterizing data in the memory 
circuit 114. In one example, such battery capacity characterizing data is stored in a 
lookup table. In another example, such battery capacity characterizing data is stored 
as a fitted or other equation. In either case, the adjusted voltage, Vadj, input 
parameter yields a resulting battery capacity reading of how much battery capacity 
has been used (or how much battery capacity remains). At 210, the battery capacity 
reading is output to either indicate how much battery capacity has already been 
used, or alternatively, how much battery capacity remains unused in the battery 102. 

FIG. 3 is a flow chart, similar in certain respects to the flow chart of FIG. 2, 
illustrating generally, by way of example, but not by way of limitation, an 
alternative technique for obtaining a reading of how much battery capacity has been 
used (or remains). At 300, the terminal voltage of the battery 102 is measured, such 
as by using the voltage measurement circuit 108. At 302, a corresponding 
measurement of load current is measured, such as by using the current measurement 
circuit 110. At 304, an actual load impedance value is calculated, such as by 
dividing the measured battery terminal voltage by the measured load current using 
the processor 112. At 306, the previously stored characterizing capacity vs. voltage 
battery data (obtained from a separate characterizing battery that is substantially 
identical in chemistry and manufacture to the battery 102) is translated to make it 
comparable to the measured battery terminal voltage, which is then used as an input 
parameter into the translated battery data. 

In one example, this includes translating the voltages in the capacity vs. 
voltage characterizing data for the battery 116, including calculating a scaled 
voltage, VscU in which Vscl = Vtest + 6*log (R1/R2), for each voltage, VtesU in the 
battery 116 characterization capacity vs. voltage, wherein R2 is a value of the 
characterizing load impedance 118, and Rl is a value of the actual load impedance 
104, and b is a coefficient obtained from characterizing the battery 116. 

At 308, the translated capacity vs. voltage data is read using the battery 
voltage measured at 300 as an input parameter. This yields a corresponding battery 
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capacity reading of how much battery capacity has been used (or how much battery 
capacity remains). At 310, the battery capacity reading is output to either indicate 
how much battery capacity has already been used, or alternatively, how much 
battery capacity remains unused in the battery 102. 

Experimental Data Analysis and Calculation of Coefficient "b " 
The above battery capacity determination techniques were analyzed for a 
pacemaker battery using lithium as the anode active material and carbon 
monofluoride (CFx) as the cathode active material. This CFx battery offers similar 
volumetric energy as lithium/iodine batteries. The CFx battery also includes, 
among other benefits, the capability to operate at medium current rates (e.g., 
between about 50 microamperes and about 500 microamperes) without suffering 
any substantial loss of battery capacity. Data was analyzed from 96 individual 
WGL Model 9424 CFx cells tested using a constant load resistance discharge. A 
promising state-of-charge model was found. The voltage-based model permits an 
implantable CRM device to provide an elective replacement indicator ("ERI," 
indicating when the implantable CRM device should be replaced due to a depleted 
battery), such as over the range from 50% to 90% state of charge of the battery. 

First, the inventors analyzed test data that included 12 cells each of WGL 
Model 9424 CFx batteries, which were placed on nine different constant current 
loads, for a total test population of 96 batteries. FIG. 4 is a graph of average battery 
terminal voltage (in millivolts) vs. battery capacity (in milliampere-hours) for each 
of the nine different substantially constant current loads. A typical DDD pacemaker 
draws an average current of about 20 microamperes, which, for the battery terminal 
voltages involved, can be approximately represented as an effective 140 kQ load 
impedance. Similarly, other load currents can be represented as effective load 
impedances (using Ohm's law) for the battery terminal voltages involved. One task 
included looking for relationships between operating voltage ("battery terminal 
voltage") and various loads. 
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FIG. 5 is a graph of battery terminal voltage (millivolts) vs. load resistance 
(kilo-ohms), as derived from the data of FIG. 4 at various levels of battery capacity. 
This can be conceptualized as taking vertical slices through the voltage curves in 
FIG. 4 at successive levels of battery discharge. FIG. 5 exhibits a pronounced 
5 exponential curvature, which appears to be consistent across the successive capacity 
levels. Suspecting that an exponential function underlies the relationship, the same 
data is plotted on a logarithmic scale of load resistance in FIG. 6. 

FIG. 6 indicates that the battery terminal voltage is, at any battery discharge 
level, a substantially linear function of the logarithm of load resistance. That is, the 
1 0 slope of this linear relationship is substantially constant at levels of discharge that 
encompass most of the battery capacity. FIG. 6 suggests that the CFx cells follow 
the Tafel relationship of Equation 2: 

77 = a ± b log i (2) 

In Equation 2, rj = overvoltage (voltage difference from open circuit) and i = load 
current (which can be represented by an approximate load impedance). Equation 2 
15 states that overpotential is a logarithmic function of load current (or load 
impedance). 

Equation 2 was used as a fitting function for the voltage-load data at several 
depths of discharge. FIG. 7 is a battery terminal voltage (millivolts) vs. load 
impedance (kilo-ohms) graph of four representative fitted curves using Equation 2. 

20 The statistical fit parameters for these linear regressions are very strong. The pure 
logarithmic behavior of operating battery terminal voltage is reflected in the very 
high degree of fit. Thus, Equation 2 is a good model for the CFx battery. The 
slopes of the fitted curves are also listed in FIG. 7. This slope is the millivolt change 
to apply in adjusting an observed operating battery terminal voltage at one load 

25 resistance to the battery terminal voltage at a second load resistance level. The 

parallelism of the lines in FIG. 7 indicates that the physical process behind Equation 
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2 is similar at each of these levels of discharge. The slope appears substantially 
constant for the continuous levels of discharge between at least 100 milliampere- 
hours and 1000 milliampere-hours. The precise linear fit and constant slope implies 
that the relationship of Equation 2 is a good physical model that applies to a wide 
5 range of loads over the major portion of the battery's service life. This allows the 
voltages reconciled for different loads to be compared over the life of the battery. 

Continuing with the regressions, Table 1 lists the fit coefficients and 
statistical parameters for several levels of discharge encompassing the entire 
capacity range of the cell. 

10 

Table 1: Fit Parameters for V=a + b log (load impedance in kilo-ohms), at various 
depths of discharge (DOD). 



standard error standard error 
coefficient (a) coefficient (b) R-square model coefficient (b) 



DOD 


Capacity 


(a) 


(b) 




s 


se(b) 


0% 


0 


2683.5 


197.69 


99.0% 


12.4 


5.86 


4% 


50 


2741.0 


110.04 


100.0% 


1.13 


0.670 


8% 


100 


2747.7 


110.18 


100.0% 


1.13 


0.537 


17% 


200 


2754.1 ! 


109.94 


100.0% 


1.13 


0.535 


25% 


300 


2752.3 


110.45 


100.0% 


1.65 


0.782 


33% 


400 


2748.2 


109.05 


99.9% 


2.21 


1.200 


50% 


600 


2752.1 


108.71 


99.9% 


2.3 


1.490 


67% 


800 


2687.8 


106.49 


100.0% 


0.17 


0.544 


83% 


1000 


2611.6 


109.52 


99.9% 


1.84 


2.820 


88% 


1050 


2580.2 


108.65 


100.0% 


0.67 


1.490 


92% 


1100 


2500.3 


126.04 


99.9% 


1.54 


3.430 


100% 


1200 


1903.7 


163.70 


91.0% 


23.8 


36.5 



1 5 The relationship of Equation 2 is to be used to adjust the measured battery 

terminal voltage for a battery 102, at some given actual load impedance 104 value, 
to that of a standard characterization load impedance 118 value for a separate state- 
of-charge characterization battery 116. Thus, we are interested in the slope 
coefficient "6" and its associated standard error. In this example, the slope 

20 coefficient "6" is seen to be consistent over a range of capacity, such as from about 
50 milliampere-hours to about 1050 milliampere-hours. Therefore, the average 
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slope, b = 109.23 mV, over that entire capacity range can be used to accurately 
translate battery terminal voltage values, at different load impedances, using a single 
function. The standard error is estimated to be 1.12 mV by averaging the standard 
errors of the slope coefficients. The calculation of voltage adjusted from resistance 
R } to R 2 including a 95% confidence interval becomes: 



In Equation 3, V adj and V obs are the adjusted (i.e., translated measured battery 
voltage) and observed voltages respectively. In this example, this single function of 
Equation 3 is valid from 4% to 88% depth-of-discharge and from 4.5 kQ to 140 kQ. 
Data for the high resistance loads is still to be collected, but is expected to be 
consistent with the data analyzed to date. 

We can test the model Equation 3 by predicting the average voltage data 
from Figure 4. A predicted voltage curve for the 32.4 kQ load was generated using 
the 16.5 kQ actual data. This is done simply by taking every point in the 16.5 kQ 
data curve and adding the term: 



This produces a predicted voltage curve for the 32.4 kQ load. Plotting predicted 
values of voltage along with the actual values provides a view of the accuracy of the 
prediction. 

FIG. 8 is a graph of battery terminal voltage (millivolts) vs. capacity 
(milliampere-hours) using both predicted and actual data for a 32.4 kQ load 
impedance. The accuracy of the fitting function is evident in the close fit of the 
projected 32.4kQ data with the actual 32.4kQ data. In this example, the fit is closest 
in the range from 50 to 1 100 milliampere-hours. Precision of the predicted voltage 
curves was evaluated by examining the residual values between the predicted and 
actual data. 
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V adj =V ob5 +[109.23 ±2.2]log 




(3) 




(4) 



FIG. 9 is a graph of the residuals (millivolts) vs. capacity (milliampere- 
hours) for the 32.4 kQ predicted and actual data of FIG. 8. The residuals for the 
32.4 kQ data reflect the precision of the fitting function. The zero line represents a 
perfect fit of the data. In the capacity range from 50 to 1 100 milliampere-hours, the 
average residual is -0.33 mV and the standard error is 1.7 mV. The range of the 
residuals is ± 3.3 mV. 

FIG. 10 is a graph of battery terminal voltage (millivolts) vs. capacity 
(milliampere-hours) for predicted and actual data for several higher resistance loads. 
The fit appears to be very good at each load. The entire set of predicted voltage 
curves were generated using the single linear operator demonstrated in Equation 4. 
FIG. 10 also provides predicted voltage curves for data that is yet to be obtained at 
the high resistance loads. In the absence of a significant self-discharge effect, these 
predicted high resistance load curves should be very accurate. The self-discharge 
rate of commercial CFx cells is believed to be very low, at about 0.5% per year. 
Therefore, the expected effect of self-discharge is practically negligible. 

FIG. 11 is a graph of the residuals (millivolts) vs. capacity (milliampere- 
hours) for the to 100 kQ predicted data (which was predicted using both 32.4 kQ 
actual data and 64.9kQ actual data) and actual data. The graph of FIG. 11 shows 
that there is no significant difference in choosing the higher 64.9kQ data as the 
predictor. Both sets of data provide an equally good fit to the actual lOOkQ values. 
This is consistent with the high degree of linearity of the slope b. 

To assess the precision of predicting battery capacity (e.g., used or 
remaining) from a discharge voltage curve, nonlinear regression was applied to the 
32.4 kQ discharge data (which was chosen as the standard because it was the 
highest load for which a completed discharge curve was obtained). The nonlinear 
regression was applied to the region extending from 600 to 1 100 milliampere-hours. 
For this 32.4 kQ load, all twelve cells were used in order to maximize the statistical 
strength of the regression. Capacity was used as the output variable since this is the 
quantity we want to predict. 
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FIG. 12 is a graph of capacity (milliampere-hours) vs. battery terminal 
voltage (volts), illustrating both the resulting cubic regression line and the lines 
indicating the 95% prediction interval. The centerline in FIG. 12 is the regression 
line and the upper and lower curves are the 95% prediction interval upper and lower 
5 bounds. The fit is very strong with an R 2 of 99.8%. All of the regression terms test 
significant. Therefore, the cubic equation is appropriate for statistically modeling 
this portion of the discharge curve, such as the example cubic equation of Equation 
5: 

Capacity (mAh) - 732235 - 814455*V + 302761*V 2 - 37557.1*V 3 (5) 

In Equation 5, V is the voltage obtained using the test load impedance 118, which in 

10 this example had a value of 32.4 kQ . The variation of voltage from cell to cell is 
small enough to allow good discrimination of battery capacity. Roughly estimating 
the average prediction interval, for a particular voltage value, the capacity can be 
estimated to within ±10 milliampere-hours. 

In view of the above experimental data using CFx batteries, the following 

1 5 example describes how to obtain a battery capacity reading. First, a unique battery 
terminal voltage reading is obtainable if the battery is discharged beyond 50% of its 
capacity. For the CFx batteries used, this corresponds to the battery terminal 
voltage dropping below about 2.95 Volts. Next, a stable battery terminal voltage 
reading is obtained, along with a corresponding stable load current reading (e.g., 

20 either directly, by averaging, or by any other technique). Then, load resistance is 
calculated, using Ohm's law, from the measured battery terminal voltage and the 
measured load current. Next, a voltage translation is performed to calculate the 
battery terminal voltage adjusted to the 32.4 kQ load impedance for which battery 
capacity vs. voltage data was previously obtained from a separate characterizing 

25 battery 116. An example of this voltage translation is indicated in Equation 6, for a 
cell which is running at 2.850 Volts and load current of 20 jxA, which is equivalent 
to a load impedance of about 142.5 kQ. 
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V adJ = 2850 + 109 log 



32.4 
142.5, 



I = 2780 mV (6) 



Using the translated battery terminal voltage, Vadj = 2.78 V, obtained from 
Equation 6, as an input parameter into the regression-derived capacity vs. voltage 
relationship of Equation 5 yields a capacity of 996 milliampere-hours. For a 1200 
milliampere hour cell, this indicates that there is about 200 milliampere-hours of 

5 remaining battery capacity. When combined with an estimate of the expected 
average current drain, an estimated time-to-end-of-life is obtained. Some 
refinements in the method would include the standard errors for the observed 
voltage and current measurement, and for cell capacity, to produce an expected 
capacity range, rather than a nominal value as demonstrated above. 

10 It is to be understood that the above description is intended to be illustrative, 

and not restrictive. For example, the above-described embodiments may be used in 
combination with each other. Many other embodiments will be apparent to those of 
skill in the art upon reviewing the above description. The scope of the invention 
should, therefore, be determined with reference to the appended claims, along with 

1 5 the full scope of equivalents to which such claims are entitled. In the appended 
claims, the terms "including" and "in which" are used as the plain-English 
equivalents of the respective terms "comprising" and "wherein." Moreover, in the 
following claims, the terms "first," "second," and "third," etc. are used merely as 
labels, and are not intended to impose numerical requirements on their objects. 
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